Effects of kink and flexible hinge defects on mechanical responses of short double 

stranded DNA molecules 



Hu CheiE and Jie YarQ 

National University of Singapore, Department of Physics, 2 Science Drive 3, Singapore 117542 

(Dated: May 3, 2009) 

We predict various detectable mechanical responses to the presence of local DNA defects which are 
defined as short DNA segments exhibiting mechanical properties obviously different from the 50 nm 
persistence length based semiflexible polymer model. The defects discussed are kinks and flexible 
hinges either permanently fixed on DNA or thermally excited. Their effects on extension shift, the 
effective persistence length, the end-to-end distance distribution, and the cyclization probability are 
computed using a transfer-matrix method. Our predictions will be useful in future experimental de- 
signs to study DNA nicks or mismatch base pairs, mechanics of specific DNA sequences, and specific 
DNA-protein interaction using magnetic tweezer, fluorescence resonance energy transfer or plasmon 
resonance technique, and the traditional biochemistry cyclization probability measurements. 
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I. INTRODUCTION 

Single-molecule stretching experiments for micron size 
DNA and cyclization experiments for DNA larger than 
230 base pairs (bp) have shown that DNA can be con- 
sidered as a homogeneous semiflexible polymer with an 
averaged bending persistence length ~ 50 nm P, 0, [3, 0| • 
The bending energy is a quadratic function of local cur- 
vature in semiflexible polymer. On the other hand, the 
rigidity and local curvature of short DNA sequences are 
sensitively sequence-dependent @ . Understanding the 
micromechanics of such sequences is important since they 
regulate protein binding and DNA compaction [7|. In 
this paper, a local DNA segment is defined as a mechan- 
ical defect if its elastic property differs from the 50 nm 
persistence length based semiflexible polymer model. A 
DNA defect can arise from sequence-dependent mechan- 
ical inhomogeneity, binding of a protein, or changes in 
DNA secondary structure. A defect can be permanently 
fixed at a specific location, or subject to thermal fluctu- 
ation. Hereafter we call the former a fixed defect, and 
the latter an excited defect. In the research, we study 
the how the defect changes several mechanical behaviors 
of short DNA. The changes may be detectable by ex- 
periments, and they reflect important information of the 
micromechanical properties of the defects. 

The defects can also be introduced by protein bind- 
ing which usually creates a local bending, or modifies 
the local bending rigidity, or both [l, [{J [lCj, [H[ • Defects 
created this way are not fixed but dynamic, subject to 
thermal fluctuation of protein binding or unbinding. In 
equilibrium, the fluctuation is determined by the inter- 
action energy between the protein and DNA. In many 
cases, proteins prefer binding to some specific sequences 
with very strong interaction. In such cases the protein- 
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binding induced defects can also be considered fixed. 

Furthermore, defects may also arise from the changes 
in the secondary structures of DNA. For examples, a 
DNA nick or a mismatch bubble that contains one or a 
few mismatched base pairs may lead to a drastic change 
in the local curvature and rigidity of DNA [l2|, G3 ■ A 
DNA bubble that contains one or a few melted base pairs 
may also appear by thermal fluctuation. In fact, the de- 
fects arising from the changes in the secondary structure 
of DNA have recently been argued to affect the microme- 
chanics of sharply bent DNA greatly [3 EE EE [HI , and 
it has been argued that, in order to correctly describe the 
micromechanics of highly bent DNA conformations, the 
traditional 50 nm persistence length based semiflexible 
polymer model of DNA should be revised to include the 
excitation of such defects [IE El, EH • 

We will focus our interest in two types of defects: one 
enhances the local flexibility, and the other creates a local 
rigid bend. Hereafter we call the former a "hinge" , and 
the latter a "kink". Traditional biochemistry methods 
to study such DNA defects are mainly based on gel elec- 
trophoresis for short DNA, in which the defect-resulted 
gel shift is observed and analyzed [2(j, or based on the 
"J-factor" measurement in which the cyclization proba- 
bility of short DNA is measured [2l[ . Recently, the rapid 
development of single-molecule manipulation and imag- 
ing methods has made it possible to study the defects 
more directly. This is because the defects can introduce 
detectable changes in DNA force response, DNA shapes, 
DNA cyclization probability, and DNA end-to-end dis- 
tance distributions. In this paper, we predict the effects 
of the hinge and kink defects on these mechanical signals. 

Effects of kinks on force-extension curve of DNA have 
been discussed previously for a long DNA. It has been 
shown that there must be sufficient number of copies of 
the defects on the DNA in order to see a detectable shift 
[22l ]. It was estimated that in order to see the effects, 
the line density of the defects on DNA should be at least 
1 per 100 nm. As a result, force-extension curve is not 
suitable to study excited defects when the excitation en- 
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ergy is too large. For the same reason, to study the 
effect of a single fixed defect, the length of DNA has 
to be not longer than ~ 100 nm. Previous single-DNA 
stretching experiments were usually done for DNA of a 
few micro meters in length. For the reasons discussed 
above, although force-extension curve measurement is a 
powerful method to study non-specific interaction DNA- 
protein interactions, it is not suitable to study specific 
DNA-protein interaction when the specific protein bind- 
ing sites are rare in DNA. In general, stretching a long 
DNA cannot detect the effects of defects that are rarely 
distributed. 

In order to study specific DNA-protein interaction or 
the mechanical property of single defects using force- 
extension curve measurement, one has to work with short 
DNA and compare with theoretical prediction at short 
length scales. Although such experiments were not re- 
ported previously, they are doable in principle. A pos- 
sible design of such experiments is discussed in section 
IV. Theoretical predictions of force-extension curve were 
not reported at short lengths, perhaps due to the lack of 
experiments of stretching short DNA. Due to the impor- 
tant applications of stretching short DNA and the lack 
of corresponding theoretical predictions, we predict the 
force-extension curve of a 100 nm DNA, and the effects of 
a single fixed kink or a flexible hinge in section III. A. We 
show that the shift in extension induced by a single defect 
is sufficiently large to be detected using single-molecule 
stretching instrument. 

The presence of defects also affects the end-to-end dis- 
tance distribution of short DNA. Experimentally, several 
methods can accurately record the dynamical fluctuation 
of the end-to-end distance of DNA. The well-known fluo- 
rescence resonance energy transfer (FRET) technique is 
able to measure dynamical fluctuation of end-to-end dis- 
tance of distances shorter than 10 nm, and it has been 
used to study the property of DNA bubble containing one 
or a few mismatched base pairs (l3| . Another very use- 
ful method is the recently developed plasmon resonance 
ruler that is able to measure dynamical fluctuation of 
distances ranged between 1 — 100 nm [23[. This method 
is perhaps more powerful since it effects in a much wider 
range and it does not have the photo bleach problem. In 
order to make predictions for possible experiments using 
plasmon resonance ruler, we study the effects of defects 
on end-to-end distance distribution for a short DNA in 
section III.B. The study covers both fixed and excited 
hinge and kink, which is more complete than previously 
reported results [lj| 11, [HI] . 

As a special case of DNA end-to-end distance distribu- 
tion, DNA cyclization probability can be measured us- 
ing biochemistry methods. Experimentally, it has been 
used to study the bending and twisting rigidities of a 
sing le stranded DNA (ssDNA) gap of several base pairs 
|12| . Theoretically, the looping probability of short DNA 
has been shown sensitive to the presence of defects. The 
effects of flexible hinge excitation and fixed defects on 
looping probability of short DNA have been studied pre- 



viously [II [11, [11 13. In section III.(C-D), in addition 
to re-confirming the previous predictions, we compute 
the effects of kink excitation, the multiple defects, and 
relative positions of the defects which have not been dis- 
cussed previously. Our results show that excited kink 
has similar effect on J-factor as excited hinge previously 
studied [l8j . suggesting an alternative mechanism to ex- 
plain the unexpected j-factor of ~ 100 bp DNA reported 
in [l4| . This excited kink mechanism is different from the 
kinkable elastic polymer mechanism previously proposed 
fl9| , where the excited defect is a hinge of zero rigidity. 

In section III.D, we also predict the end meeting an- 
gle dependence of DNA looping and its interaction with 
the presence of defects. Such predictions are important 
in studies of DNA loops where the ends are bridged by 
a protein, since such DNA looping proteins may likely 
impose a specific requirement on the meeting angle of 
DNA. Previously, it has been demonstrated that the pro- 
tein size greatly affects the looping probability (26[, but 
the effects of meeting angle has not been studied. Our re- 
sults show that the looping probability of a 40 nm DNA 
is extremely sensitive to the meeting angle and to the 
presence of defects. 

The computation of the force-extension curve, the end- 
to-end distance distribution, and the cyclization proba- 
bility of DNA are done us ing the transfer-matrix methods 
previously developed in jl8l |25| . We believe our predic- 
tions will be useful for future experiments designed to 
study the micromechanical properties of the defects. 



II. THEORY AND MODEL BASED ON 
TRANSFER MATRIX METHOD 

We consider the discrete polymer model consisting of 
N segments with equal length b. In such a model, the 
conformation of DNA is described by the orientations of 
the segments, ti, where i = 1, ••• ,N, and the confor- 
mational energy is carried by the bending of the vertices 
connecting adjacent segments. The energy in unit of fc^T 
of this model is a summation of the vertex energies: 

N-l 

E = Y, E Si,U+i)- (1) 
1=1 

For semiflexible polymer model of DNA, Ei — 

^ (tj + i — tj) , where a = A/b is a dimensionless quantity 
describing the rigidity of the vertex. The parameter A is 
the persistence length of the semiflexible polymer which 
describes the local rigidity of the polymer. In order to 
approximate the continuum semiflexible model, it is re- 
quired to choose a segment length b <C A. For DNA, it 
has been accepted that A = 50 nm under normal experi- 
mental conditions [J @, 0] . Our computations described 
in the following sections are based on a choice of b = 1 
nm, thus a = 50 for DNA in the 50 nm persistence length 
based semiflexible polymer model. 
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the flexible hinge 
t, ) , where a! is 



The site-dependent vertex energy ^(tj, tj+i) is able 
to describe fixed defects. Two kinds of mechanical de- 
fects are discussed in the paper: 1 

defect Ei t hinge(ii+1, t,) = y (ti+i 

much smaller than the vertex rigidity of B-form DNA, 
and 2) the kink defect E iy kink = " ti+i — 7) 2 , where 
7 = cos((9) and is the preferred (lowest-energy) kink an- 
gle. These defect energy forms have been used in some of 
our previous studies O, H2, [H| . In the computation, we 
choose a' — 1 for a flexible hinge vertex, and (9 = tt/2, 
a=50) for a kink vertex. 

To describe the excited state, a vertex can be made 
to have two states, indicated by an index nf. ni = 
for B-form state, and rii = 1 for defect state. With the 
excitation, the vertex energy becomes: 



Ei(U, tj+i) — 5 ni fiEifi(ti, tj+i) 



(2) 



where .E^o is the B-form DNA vertex energy, and E^ is 
the defect energy. The parameter [i controls the proba- 
bility that a defect forms, and represents the free energy 
cost of generating a local defect [l8l.[25j. 

A transfer-matrix method has been developed to com- 
pute the end-to-end distance distribution, the looping 
probability, and the force-extension curve [25| . The par- 
tition function of conformations subject to certain con- 
straints can be calculated by the inverse Fourier trans- 
formation of the partition function in wave number space 
Z\l. Z\l can be computed by a direct order product of the 
site-dependent matrices, so the finite size effect of short 
DNA, sequence- dependence, and the boundary condition 
are automatically included. The probability of specific 
DNA conformation is the ratio of the constrained parti- 
tion function and total partition function. To calculate 
the extension of DNA subject to certain stretching force, 

N-l 

a force induced energy ^2 —bfz ■ tj is added into the 

i=i 

energy equation. Extension of DNA is obtained from the 
derivative of the free energy with respect to force. Please 
find the details of the transfer-matrix method in 12511. 



III. DNA MECHANICAL RESPONSES TO THE 
DEFECTS 

A. DNA extension under tensile force 

In a single molecular manipulation experiment, a ten- 
sile force can be applied to DNA and the resulted ex- 
tension of DNA can be measured. The mechanical prop- 
erties of DNA can be derived from the force-extension 
curve @,||. In such experiments, the DNA usually has 
a contour length of a few micrometers. In [22j], it was 
shown that in order to see detectable effects of defects on 
the force-extension curve, the density of the defects on 
DNA must be sufficiently large. The minimal density of 
kink defects was shown to be around one per 100 nm [22j |. 



Because of this, one does not expect to see the effects of 
excitation of defects if the excitation energy is too high 
[III], or the effects of a single fixed defect if the DNA is 
too large. In this section, we predict the detectable ef- 
fects of single fixed defects on the force-extension curve 
of a short DNA of 100 nm. 

The force-extension curve may be computed using the 
site-dependent transfer-matrix method developed in [25| . 
Fig. [1] shows the force-extension curve of a 100 nm long 
DNA that contains a single defect in the middle. From 
the figure one can see a single fixed defect results in de- 
tectable shift of the force-extension curve. Inset of Fig. 
[1] (a) shows the shift of the defect-containing DNA curve 
from the defect-free DNA curve. The largest shift is ~ 7 
nm for DNA containing a flexible hinge, and ~ 12 nm 
for DNA containing a kink, at a force ~ 0.14 pN. Such 
shifts are detectable in single molecule experiments. Al- 
ternatively, one can study the effects of the defects by 
measuring the effective persistence of DNA. The effec- 
tive persistence length may be determined by fitting the 
well-known force-extension curve at large force approxi- 
mation 0: < z > /L = 1- ^k B T/4Af, where < z > / L 
is the average extension scaled by the contour length, / 
is the tension along z direction, and A is the effective 
persistence length. Fig. QJb) shows the results: effec- 
tive persistence of the defect-free DNA determined by 
the slope fitting is A ~ 36 nm, apparently smaller than 
the expected 50 nm. This reduced value is due to the 
finite-size effect which was recently discussed in [27| • The 
slope fitting for defect-containing DNA gives A ~ 30 nm 
for DNA containing a single flexible hinge, and A ~ 22 
nm for DNA containing a single 90° kink. Longer DNA 
with contour length 500 nm is also studied, where the 
effects of the defects become much less detectable. The 
effective persistence length is A ~ 46 nm for defect-free 
DNA, A ~ 44 nm for DNA containing a flexible hinge, 
and A ~ 41 nm for DNA containing a 90° kink. 

When tension approaches zero, the average extension 
of DNA also approaches zero. As such, at small tension, 
the radius of end-to-end distance (r) is a better measur- 
able quantity [28| . Subject to a tension along z direction, 
< r 2 > is determined by the following equation: 



< r 2 >=< z > 2 +aj + - 



f 



(3) 



where a 2 is the variance of DNA extension which can 
be obtained from the secondary derivative of free energy. 
The effects of defects on < r 2 > are shown in Fig. O At 
small force, < r 2 > does not vanish as expected, and it is 
determined by the radial end-to-end distance distribution 
of a tension-free DNA which will be discussed in the next 
section. 



B. End-to-End distance distribution 

An approximate analytical expression has been derived 
for the radial end-to-end distance distribution of short 
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FIG. 1: The effects of fixed defects on force-extension curve 
measurements, (a) The force-extension curve of a 100 nm 
DNA that does not contain defects (solid curve), contains 
a single flexible hinge defect at center (dashed curve), and 
contains a single 90° kink defect at center (dotted curve). 
Inset shows the magnitude of the shift of the defect-containing 
DNA curves from the defect-free DNA curve, (b) The slopes 
are used to compute effective persistence lengths, which are 
found to be ~ 36 nm for the defect-free DNA, ~ 30 nm for 
the fixed flexible hinge containing DNA, and ~ 22 nm for the 
fixed kink containing DNA. 



FIG. 2: < r > as functions of stretching force for defect-free 
and defect-containing DNA. The symbols are the same as in 

Fig. m 



the excitation of the defects are much less obvious in Fig. 
0(a). This is because the distribution is most sensitive to 
the defects only when R is very small where the excited 
defects begin to dominate the bending. To see it, we 
rcplot the distribution using logarithm scale as shown in 
Fig. [3] (b). Our results show that, when R < 10 nm, 
47ri? 2 p(R) for the DNA subject to excitations of defects 
is about two orders of magnitude larger than that of a 
defect-free DNA. 

The position of the fixed defect affects the end-to-end 
distance distribution. We computed the positional effect 
of a fixed 90° kink and a fixed flexible hinge for a 100 nm 
long DNA. For both defects, their effects are maximized 
when they are in the middle of DNA. When the defect 
moves from one end to the center of the DNA, the peak 
of the distribution moves to shorter end-to-end distance, 
and the distribution becomes wider. This observation is 
in agreement with the computation done in [24] where 
the positional dependence of a fixed flexible hinge was 
studied. 



scmiflexible polymers [29| . By using the transfer matrix 
method developed in [25[ , the distribution of the end-to- 
end distance vector (R) of a polymer subject to defect 
excitation or containing fixed defects can be calculated. 
In many cases the vector distribution p(R.) does not de- 
pend on the direction of R, and in such cases the scalar 
end-to-end distance distribution is simply AirR 2 p(R). 

Fig. [3] shows how fixed and excited defects affect the 
scalar end-to-end distance of a 40 nm DNA. In the studies 
of the excitation of defects, the excitation energy for kink 
is 10 fcgT, and that for flexible hinge is 11 fceT. Clearly, 
the distribution of end-to-end distance of short DNA is 
extremely sensitive to the fixed defects as shown in Fig. [3] 
(a) . The peak of the distribution moves to smaller end-to- 
end distance, and becomes wider. However, the effects of 



C. J-factor measurement of DNA molecules 

The previous section shows the effects of defects on 
the scalar end-to-end distance distribution 47ri? 2 p(R). 
Their effects are found much enhanced when R is small, 
and one may expect the effects are maximized at R = 
where DNA loop forms. Experiments can be carried out 
in such a way that the rate of cyclization of DNA is pro- 
portional to p//(0) which is just p(0) under a constraint 
that the two DNA ends are parallel to each other [30l |. 
The quantity measured in such experiments is the "J- 
factor", which is related to p//{0) by J = jfcp//(0), 
where Na is Avogadro's number. This quantity, usually 
expressed in units of mol/liter (M), is a measure of the 
equilibrium concentration of one end of the DNA at the 
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FIG. 3: Effects of defects on scalar end-to-end distance dis- 
tribution of a 40 nm long DNA plotted in linear scale (a) 
and logarithm scale (b). In (a), curves for excited defects are 
nearly overlapped with defect-free DNA. The fixed defects are 
located in the middle of the DNA. For excited defects, exci- 
tation energy for kink is 10 ksT, and that for flexible hinge 
is 11 k B T. 



position of other end. Several methods have been devel- 
oped to compute the J-factor for semiflexible polymers 

In the transfer-matrix computation, due to the parallel 
boundary condition, the constrained partition function of 
circular conformations is the trace of the ordered product 
of the site-dependent matrices [25[ . Because of the cyclic 
property of trace, the cyclization probability of DNA con- 
taining a single fixed defect does not depend on the lo- 
cation of the defect. The details of the transfer-matrix 
computation of DNA cyclization probability is described 
in [23. The effects of excitation of flexible hinge and a 
single fixed kink on J-factor are discussed there. In this 
section we report a more complete computation for the 
effects of the defects, including the kink excitation and 
multiple defects which have not been discussed earlier. 

Fig. E] shows the effects of fixed or excited defects 
described previously on DNA J-factor. For long DNA, 
the cyclization of DNA is dominated entropically, so the 



FIG. 4: The effects of the positions of fixed defects on the 
scalar distribution of end-to-end distance of a 100 nm long 
DNA. (a) The positional effect of a fixed 90° kink, (b) The 
positional effect of a flexible hinge. In (a) and (b), the ref- 
erence distributions of the defect-free DNA are displayed in 
gray color. The dark curves, from top to bottom, are results 
computed for defect locations at 5, 10, 20, 30, 50 nm from 
one end of DNA, respectively. 



rarely distributed defects do not affect the J-factor much. 
On the other hand when DNA is shorter than ~ 100 
nm, DNA bending energy begins to dominate the cycliza- 
tion of DNA, and the effects of defects becomes enhanced 
since they greatly reduce the free energy of DNA circles. 
From Fig. [5j the effect of a fixed defect becomes notice- 
able when DNA is less than around 200 nm, and becomes 
greater when DNA is shorter. Comparing with a single 
fixed defect, the effect of an excited defect becomes no- 
ticeable at a shorter length (< 80 nm). This is because 
there is a penalty of excitation energy of 10 ksT for kink 
or 11 ksT for hinge to have an excited defect. When 
DNA is shorter than ~ 40 nm, DNA with excited defects 
even begin to have J-factor larger than DNA with a sin- 
gle fixed defect, indicating that more than one defect is 
excited to further relax sharp bending of the mini circular 
DNA. Interestingly, J-factor of DNA with excited flexi- 
ble hinge defect increases with decreasing DNA contour 
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length when it is shorter than 40 nm. This phenomena 
can be easily understood: when there are two flexible 
hinges already excited at optimal positions to allow a 
direct folding to form DNA loop satisfying the parallel 
boundary condition, the looping is dominated by how 
easy the two nearly "freely-joint" segments can find each 
other. Obviously, the smaller the circle size, the larger 
the chance the two ends meet. As shown in [ill Ojj], ex- 
cited flexible hinge can give high J-factor for short DNA 
circles. Fig. [5] shows that excited rigid 90° kink also has 
similar effect. 

When there are more than one fixed defects along 
DNA, J-factor of such a polymer then depends on the 
locations of the two defects. Fig. [5] shows the how J- 
factor depends on the separation between two fixed flex- 
ible hinge defects on a 50 nm DNA. Clearly, J-factor is 
very sensitive to the relative position of the two defects. 
It increases about 4 orders of magnitude when the sep- 
aration between two defects increases from near zero to 
half of the circle size, 25 nm. The optimal separation 
is half of the circle size, because this separation allows 
a direct folding for loop formation while satisfying the 
parallel boundary condition. The looping probability de- 
pends only on the separation between the two defects in 
the loop conformation, which is also a result of the cyclic 
property of trace. 
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FIG. 5: J-factors as functions of DNA length. The fixed 
defects are located in the middle of DNA. For excited defects, 
the excitation energy for kink is 10 fcsT, and that for flexible 
hinge is 11 fcsT. 



D. Looping probability under free boundary 
condition and dependence on meeting angle 

The J-factor measurement discussed in the previous 
section essentially measures the cyclization probability 
with a parallel boundary condition. DNA loops can also 
form under other boundary conditions. For example, a 
protein mediated loop may impose a specific meeting an- 
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FIG. 6: The effect of separation of two fixed flexible hinge 
defects on J-factor of a 50 nm DNA. 



gle constraint on the two DNA ends. If the interaction of 
the ends does not impose a strong orientational require- 
ment, then we call it "free boundary condition" . The 
effects of excitation of flexible hinge and a single middle- 
located fixed kink under the free boundary condition were 
also discussed in [25| . The effect of the position of a sin- 
gle fixed flexible hinge on looping probability of a short 
DNA under free boundary condition was computed in 
[24l ]. Here we report a more complete computation in- 
cluding the effects of fixed flexible hinge on a larger length 
range, the excitation of kinks, the positional dependence 
of kink, and the dependence on meeting angle constraint 
which have not been discussed earlier. 

As shown in Fig. an interesting observation is that 
the looping probability of DNA with a middle-located 
fixed flexible hinge scales with loop length with a power 
law. This is because the hinge basically breaks the short 
DNA into two nearly freely joint segments. If the hinge 
is completely free to rotate and the two identical arms 
are infinite rigid, the looping probability scales inversely 
with the volume which the ends can sample. As a result, 
one expects Pi 00 ping °c 1/L 3 . In our computation, the 
flexible hinge has a finite bending rigidity and the two 
arms are not infinite rigid, leading to a slightly different 
scaling: Pi oop ing °c 1/L 3 5 . The curve corresponding to 
the excitation of flexible hinges displays similar power law 
when DNA length is shorter than 40 nm. This is because 
when DNA becomes very short, a flexible hinge is always 
excited in the middle of the DNA to allow a direct folding 
to meet the two ends. The curve of the fixed kink does 
not display a simple power law. The large rigidity of 
the kink (a = 50) leads to decreasing looping probability 
when the length is getting shorter than ~ 40 nm. 

Fig. [5] shows the looping probabilities as functions of 
the meeting angle for defect-free DNA, excited and fixed 
defect-containing DNA with a contour length of 40 nm. 
For a perfect semiflexible polymer, the looped conforma- 
tion with lowest bending energy is like a teardrop with 
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meeting angle around 1.8 rad [34|. This is confirmed 
by our calculation: the looping probability of defect- 
free DNA is maximal when meeting angle is around 
1.8 rad. The defects increase the looping probability 
under anti-parallel boundary condition (Cos(9) = — 1) 
the most while that under parallel boundary condition 
(Cos(8) = 1) nearly the least. This is because a single de- 
fect permanently located or excited at the middle is suf- 
ficient to allow DNA loop formation under anti-parallel 
condition, while under the parallel condition, it requires 
up to 2 defects to allow DNA to be looped. As such, 
under free boundary condition, defected DNA tends to 
form loops with anti-parallel meeting angle. In the pres- 
ence of non-specific DNA crosslinking factors, such defect 
formation favors the formation of hairpin structures. 

Unlike the parallel boundary condition, looping under 
free boundary condition depends on the location of a sin- 
gle fixed defect. This is because the partition function is 
not the trace of the product of the matrices under the free 
boundary condition, as shown in 25]. Fig. [5] shows the 
effects of the locations of single kink and single flexible 
hinge defects for looping of a 100 nm DNA. In agreement 
with the observation in [24], the optimal defect position 
is the center of DNA. 




Loop length (nm) 




FIG. 8: The dependence of looping probability on the meet- 
ing angle between the first and last tangent vector for 40 
nm DNA. Cos(8) — 1 for parallel boundary condition, while 
Cos(8) — —1 for anti-parallel boundary condition. The sym- 
bols are the same as in Fig. [7] 
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Position of defect along 100 nm DNA (nm) 



FIG. 9: Dependence of looping probability of a 100 nm DNA 
on the positions of single fixed defects. 



FIG. 7: Effects of defects on DNA looping under free bound- 
ary condition. The fixed defects are in the middle of the DNA. 
For excitation of defects, the excitation energy for kink is 10 
UbT and that for flexible hinge is 11 UbT. 



IV. DISCUSSION AND CONCLUSION 

We have computed the effects of kink and flexible hinge 
defects on various mechanical responses of short DNA. 
We showed that the defects can lead to detectable shift in 
extension of a ~ 100 nm long DNA subject to a constant 
tension using a magnetic tweezer. We also showed that 
the effects of fixed defects can be detected by the changes 
in effective persistence lengths. To obtain the persistence 
length, one has to obtain accurate measurement for the 



tension and the extension. In magnetic tweezer, the ten- 
sion is determined by the transverse fluctuation of one 
end of DNA and the extension of DNA. For short DNA, 
the relative errors in the measurements of both quantities 
usually become larger, leading to some potential techni- 
cal difficulties to measure the effective persistence length 
with high accuracy. One possible method to overcome 
this difficulty is to use a long DNA handle which tethered 
between a paramagnetic particle and a non-paramagnetic 
particle. The non-paramagnetic particle is attached to 
the end of the short DNA being studied. The force is 
now applied to the paramagnetic particle, and it can be 
accurately determined using the long DNA handle. The 
extension of the short DNA is determined by the position 
of the non-paramagnetic particle. Using such a design, 
the accurate force-extension curve and the effective per- 
sistence length of a short DNA may be determined. 
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We then computed the effects of defects on the end- 
to-end distance distribution of a 40 nm DNA. Such mea- 
surements can be done by recording the dynamical plas- 
mon resonance wavelength of two 40 nm golden parti- 
cles attached to the ends of DNA [23|. This technique 
was reported able to monitor the dynamical changes in 
DNA end-to-end distance over a range of 1-100 nm. As 
such, our predictions of the end-to-end distance distri- 
bution of the 40 nm DNA are experimental detectable. 
As shown in Fig. [3J both the fixed and excited defects 
affect the distribution, particularly when the end-to-end 
distance is shorter than 10 nm. If the DNA being studied 
is shorter than 10 nm, then the end-to-end distance fluc- 
tuation can also be measured by FRET technique. The 
dynamical end-to-end distance fluctuation measured by 
plasmon resonance or FRET can be used to study the dy- 
namical binding of protein to DNA, the mechanical prop- 
erty of specific DNA sequences (for example, the nucleo- 
some positioning sequences) , the mechanical property of 
highly bent DNA, and other DNA defects such as DNA 
nicks and DNA mismatches. 

Although there have been many theoretical computa- 
tions for effects of defects on DNA cyclization and looping 
probability, we still include a section discussing it in or- 
der to make the research more complete. In addition to 
confirming a few previously reported predictions, some of 
our computations have not been discussed before. These 
include the effect of a fixed flexible hinge over a larger 
length range, the effect of multiple defects and excita- 
tion of rigid kinks, as well as the dependence on meeting 
angle. In the computation of the effect of a fixed defect 
on DNA J-factor measurement, one can see in Fig. [5] 
that a single fixed flexible hinge takes a detectable effect 
up to loop size of ~ 200 nm (600 bp). We note that 
the effect of a single fixed flexible hinge on J-factor mea- 
surement is equivalent to the effect of a flexible meeting 
angle boundary condition that imposes an energy of the 
meeting angle (9) of the ends, E{9) — ^9 2 , on J-factor 
measurement of a defect-free DNA of the same length. 
Thus, one would predict that the effect of such an elastic 
meeting angle constraint can affect the J-factor up to a 
looping size ~ 200 nm if taking a = 1. Recently, the flex- 
ible meeting angle boundary condition was proposed [35[ 
to explain the discrepancy between the J-factor measure- 
ments of ~ 100 bp DNA reported in [l4[ and (36|. In this 
mechanism, the possible effect of the two nicks after the 
annealing of two DNA overhangs and before the ligation 
reaction is considered. The nicked, annealed overhangs 
of 4 bp were assumed to be a flexible hinge in the mech- 
anism. By varying the rigidity of the defect, it is shown 
that this mechanism is able to explain the both experi- 
ments using a = 1 for J-factor reported in [14| and a — 20 
for J-factor reported in [l4|]. Noticing the equivalence 
between this mechanism and the fixed flexible hinge in 
their effects on J-factor, therefore, this mechanism does 
not explain J-factor measurements for DNA larger than 



230 bp, where all the measurements are in good agree- 
ment with predictions by 50 nm persistence length based 
semiflexible model [lElli]. 

The effects of meeting angle constraint on DNA loop- 
ing have not been discussed previously. In Fig. [5J we 
show that the looping probability of a 40 nm defect-free 
DNA is sensitive to meeting angle, with a variation of 
about two orders of magnitude depending on the value of 
meeting angle. The sensitivity of DNA looping to meet- 
ing angle has important biological implications, since it 
is required to understand the looping efficiency of DNA 
with the ends bridged by proteins that impose a specific 
constraint on the meeting angle. As mentioned in intro- 
duction, another important factor affecting DNA looping 
efficiency is the size of the protein bridge [26| . We also 
show that, in the presence of defects, the looping prob- 
ability in the whole angle range can be shifted up by at 
least two orders of magnitude at 9 = 0, and up to seven 
orders of magnitude at 9 = tt. This suggests that DNA 
with a defect tends to form loops in anti-parallel config- 
urations. It implies formation of hairpin structure when 
the loop is stabilized by DNA bridging proteins such as 
HN-S. A defect would initiate formation of a seeding loop 
with a near anti-parallel boundary condition which is sta- 
bilized by a single bridging protein, followed by progress- 
ing into a hairpin structure by adding more proteins se- 
quentially. 

All our computations described in the paper are done 
based on the transfer-matrix method developed previ- 
ously [25[ • This method allows us to study the sequence- 
dependent properties, fixed or excited defects, and dif- 
ferent boundary conditions. There is, however, a draw- 
back in the current transfer-matrix method: there is no 
easy way to include the twist constraint into this method, 
since with twist we no longer have localized interaction 
which only depends on adjacent orientations of DNA. 
Since twist constraint affects DNA looping probability 
and shapes of DNA loops greatly for small DNA, one of 
our future researches is to find a way to include the twist 
constraint into the transfer-matrix computation. 

Finally, we believe our predictions will find its applica- 
tions in designs of future experiments to study DNA de- 
fects, specific DNA sequences, and specific DNA-protein 
interaction, using magnetic tweezer, FRET or plasmon 
resonance techniques, and the traditional biochemistry 
cyclization probability measurement. 
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